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Ph. Abstract 

CD . 

'-^ ' Recently it has been realized that the production and decay processes of charginos, neu- 

tralinos, and sleptons receive corrections which grow like logrriq for large rriq. In this paper 
we calculate the chargino pair production cross section at e + e~ colliders with quark/squark 
loop corrections. We introduce a novel formulation, where the one-loop amplitude is reor- 
ganized into two parts. One part is expressed in terms of the "effective" chargino coupling 
g e ~^ and mixing matrices U p , V p , and includes all 0(logm,q) corrections, while the other 
decouples for large mg. The form of the one-loop cross section then becomes physically trans- 
parent. Our formulation can be easily extended to other loops and processes. Numerically, 
we find significant corrections due to the effective t-channel coupling g e -^, for gaugino-like 
charginos. In the mixed region, where the chargino has large gaugino and Higgsino compo- 
nents, the corrections due to (U p , V p ) are also significant. Our numerical results disagree 
with a previous calculation. We revisit previous studies of the determination of 5 e -^> through 
the measurement of the chargino production cross section. We point out that a previous 
study, which claimed that the measurement suffers large systematic errors, was performed 
at a "pessimistic" point in MSSM parameter space. We provide reasons why the systematic 
errors are not a limiting factor for generic parameter choices. 



1 Introduction 



Calculations of higher loop effects in the Standard Model (SM), together with the recent precision 
measurements of electroweak parameters, have given rise to a wealth of information on physics at 
the weak scale and above. Among these measurements, one of the interesting observations is the 
approximate agreement of the measured gauge couplings with the prediction of supersymmetric 
Grand Unified Theories (GUTs) Jl|]. This may be regarded as indirect evidence for the minimal 
supersymmetric standard model (MSSM), which is the low energy effective theory of supersym- 
metric GUTs. Additionally, global fits to precision data in the SM prefer a light Higgs boson mass 
||, which is consistent with the MSSM which predicts < 130 GeV. 

In the near future, the Large Hadron Collider (LHC) will explore the TeV energy region. 
Squarks and gluinos will be discovered, together with charginos and neutralinos, if the supersym- 
metry breaking scale is below a few TeV. Recent studies show that certain superpartner mass 
differences can be measured quite precisely at the LHC || ||. Furthermore, if any of the proposed 
lepton colliders are constructed many of MSSM parameters (e.g. the gaugino masses, the Higgsino 
mass /x, the slepton masses, and the ratio of vacuum expectation values tan/3) will be measured 
to 0(1%) ~ 0(10%) |, |, 0,1- 

The future precision measurements of new particle masses, event rates, and branching ratios, 
will provide for detailed tests of the supersymmetry hypothesis. Supersymmetry imposes hard 
relations between gaugino couplings and gauge couplings, and between Higgsino and Higgs cou- 
plings. The cancellation of Higgs mass quadratic divergences cannot be realized without these 
supersymmetry coupling relations. Therefore, they comprise an essential ingredient of the model. 
Measurements of the coupling relations will provide definitive evidence of the realization of super- 
symmetry in nature. 

Because supersymmetry is broken, the hard coupling relations receive radiative corrections 
II, H, [10], [H]? [12], [I3l| . Since all split supersymmetry multiplets contribute to the splitting of the 



gauge/gaugino and Higgs/Higgsino couplings, measurements of the splitting may provide useful 
information about the supersymmetry spectrum. This is readily understood from the point of 
view of effective field theory. As an example, below the squark mass scale the gauge and gaugino 
couplings run differently because squarks do not contribute to the running of the gauge or gaugino 
couplings, but quarks continue to contribute to the running of the gauge couplings. At a scale 
Q below the squark mass scale, this mismatch in the running manifests as a difference between 
the couplings proportional to \n{rriq/Q) ||. Such a correction also appears in the off-diagonal 
elements of the chargino and neutralino mass matrices, which originate from Higgsino-Higgs- 
gaugino couplings. Notice the analogy to the radiative corrections in the SM. The SM gauge 
symmetry relates various SM observables. Since the particle masses in the SM break the gauge 
symmetry, the measurement of electroweak observables leads to constraints on the particle masses, 
in particular the top quark and Higgs boson masses. 

Corrections to supersymmetric relations were first calculated in Ref. || in the effective renor- 
malization group equation (RGE) approach. For degenerate squarks, one finds the correction to 
lepton-slepton-gaugino couplings 

9 e ,wl9l M = 1 + 2% log 10 K~/m r ) (1) 
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g^/gT = 1 + 0.7% log 10 K~/m f ) . (2) 

The fermion-sfermion-gaugino couplings are involved in both the production and decay pro- 
cesses of charginos, neutralinos and sfermions. Studying these processes provides for measurements 
of the gaugino couplings. It is particularly interesting to measure the gaugino coupling at future 
e + e~ colliders, because precise measurements of the differential cross section are possible there. 
Studies show the gauge/gaugino coupling difference may be measured within 0.3% ~ 20% through 



the measurement of the production cross sections of sleptons or charginos || [11], [L3|]. Typically, 
the high sensitivity of (9(1%) or less may be achieved when the collider experiments can measure 
both the final state superpartner masses and the mass of the particle exchanged in t-channel. Such 
a high precision measurement of the difference between the gauge and gaugino couplings allows 
for the possibility of constraining the mass scale of squarks which might not be in direct reach in 
either hadron or lepton collider experiments. 

Because the corrections to the supersymmetric relations are large enough to be measured in 
proposed future experiments, it is important to calculate the full one-loop amplitude in detail. 



Tree level amplitudes depend on the definition of the tree level parameters. In the DR scheme [H 



the gauge couplings and chargino and neutralino mixing matrices depend on the renormalization 
scale. Changing the scale by a factor of 2 easily results a few percent change in the predicted value 
of the production cross section, about the size of the correction of interest. Such scale dependence 
can be curtailed only by including radiative corrections. 

In this paper we present the full one-loop calculation of the chargino production cross sec- 
tion cr(e~e + — > XiXj) including quark and squark loop contributions. The calculation has been 
performed previously in the DR scheme in Ref. |D|. In their formula the mixing matrix of the 



chargino is scale dependent. This scale dependence must be compensated for by the chargino wave 
function renormalization, leading to very complicated expressions. We simplify the calculation by 
introducing the effective mixing matrices U p , V p . Expressed in terms of U p , V p , the formulae are 
reorganized into a compact and physically transparent form. This reorganization allows us to see 
that the full amplitude consists of two renormalization scale independent parts. One contains all 
the process independent corrections. For sufficiently heavy squarks, this reduces to the effective 
tree level amplitude which depends on process independent effective parameters. These effective 
parameters contain all the corrections proportional to logrriq. The other part of the amplitude 
contains the process dependent contributions, i.e. the one particle irreducible (1PI) chargino ver- 
tex correction and chargino wave function renormalization. This part decouples in the large rriq 
limit. 

In this paper, we also examine previous studies of the measurement of the effective gaugino 
coupling g~ e ~yr through the study of chargino production and decay |IT| , which is based on Monte 
Carlo (MC) study of Ref. . In the study, the constraint on g ™ is claimed to be limited by the 
systematic error due to the uncertainty of the underlying parameters. The maximal sensitivity to 
3evW obtained in Ref. is 2%, which is merely enough to constrain squark mass within a factor 
of 10. One may ask whether a full one- loop calculation is necessary if this is always the case. 
However, we find the case studied in Ref. JjJ is uncharacteristically pessimistic in the sense that 
the signature of chargino events is very similar to that of backgrounds, and this naturally makes 
precision measurement very difficult. We provide reasons why systematic errors are not a limiting 



2 



factor in the precision study of the supersymmetric relation. 

The paper is organized as follows. In Sec. 2 we describe our formalism which reorganizes 
the one-loop chargino pair production amplitude, making for a more physically transparent and 
computationally manageable formula. In Sec. 3 we show our numerical results. The formulas 
presented in Sec. 2 systematically guide the discussion of the logm^ corrections and the remaining 
finite corrections. In addition to corrections to the amplitude from the well studied logm^ behavior 
of g ei> ^Y, we find the logm^ corrections and some finite corrections to the effective mixing matrices 
are important when the chargino is a sizable mixture of Higgsino and gaugino. We also note 
in Sec. 3 that the calculation of the cross section including only top-stop and bottom-sbottom 
loops presented in Ref. |15|] is not a reasonable approximation in general. Our numerical results 



disagree with the results of Ref. ||15|| . We discuss the validity of various approximations to the 
one-loop cross section. In Sec. 4 we revisit the previous study of the g e ~^ measurement from the 



chargino production cross section. In Ref. [11] they argue that the sensitivity to the cross section 
is limited due to a strong dependence of the chargino acceptance on the theoretical underlying 
parameters. However, this study is carried out at a point in the supersymmetry parameter space 
with very special kinematics. Unlike the generic supersymmetric signature, the signal has a soft ]/>t 
distribution similar to the WW background. This results in the small acceptance and the strong 
sensitivity to the masses under the standard set of cuts. We point out that at a generic point in 
parameter space the acceptance error is not large and the systematic error is not a limiting factor 
in the measurement. We also discuss general ways to minimize the acceptance error. Sec. 5 is 
saved for discussion and conclusions. 

2 One-loop correction to chargino pair production 
2.1 Amplitude and cross section 

We show the form of the amplitudes of the chargino pair production e~(pi)e + (p2) — ► Xi {Pz)Xj (pi) 
including quark and squark loop corrections. We start with the tree level amplitude. The s-channel 
amplitude comes from the exchange of gauge bosons (7, Z). The t-channel amplitude involves the 
exchange of the electron sneutrino v = v e . Their sum gives 

iMij = i(M { ^ +Mf j > ) 

= +i e 2 - [u Os ) 7 M $ij v (p 4 ) ] [v (p 2 ) 7m m (Pi ) ] 



+^^^Hp 3 h"PRv(pdMP2h,PM Pl )} • (3) 



where s = (pi+p 2 ) 2 , t = (pi -p 3 ) 2 , 9z = W c w, s w = ^l-c 2 w = s\nO w , and P R:L = (l±7 5 )/2. 
The u and v are the wave functions of e ± and ■ We have applied a Fierz transformation to the 
t-channel amplitude to make its spinor structure similar to s-channel one. The vfj are the tree 
level couplings of the charginos to the Z boson, which depend on the chargino mixing matrices 
(U, V). Their definitions are given in Appendix A. 
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We next show the form of the corrected amplitude. The loop corrections include the 1PI 
chargino/gauge vertex correction, the chargino wave function renormalization, and the gauge boson 
self energy corrections. We adopt the DR renormalization scheme for gauge couplings and the weak 
mixing angle (e, g%, gz, sw) and chargino mixing matrices (U, V), and we adopt the on-shell scheme 
for the Z and v masses. Since we only include the quark and squark loop corrections, the running 
of the DR parameters includes only the contributions of quarks and squarks for consistency. Note 
that (U, V) are obtained by diagonalizing the tree level mass matrix Eq. ( A. 2 ) in the DR scheme, 
so they are renormalization scale dependent. 

The t-channel amplitude receives only chargino wave function renormalization. The corrected 
amplitude takes the following form 

x tysvn + hyiv^sz* + 5z*)v; 1 v jl )) , (4) 

where 5Z L ^ are the wave function renormalization of charginos Xl(r)- Their explicit forms are 
given later. Strictly speaking, the squark loop correction also appears in the sneutrino propaga- 
tor. However, this momentum independent contribution is completely canceled by the on-shell 
renormalization of m~. 

The s-channel amplitude is corrected by oblique gauge boson propagator corrections, chargino 
wave function renormalization, and the 1PI chargino-chargino-gauge boson vertex correction. The 
form of the corrected amplitude is 



iM® = -ieQe- ( 1 - ) Hp^v^Mp^^)} 




-id 1 TToT l Qx- 6 iA U (P3h fl v(p i )}[v(p 2 h^gz(hePL ~ S w Q e )u{pi)\ 



+Qe[u{pz)l»{v Z Lij P L + vf tij P R )v(p 4 )}[v(p 2 )^u( Pl )] 

• 1 L _ n T zz(s)-n T zz(M z ) \ 

9z s-Ml\ s-M\ ) 

x[u(p3)Tf j fl v(p4)][v(p2) r r^(h e PL ~ *wQe)u(pi)] , (5) 

where Mz is the Z-boson pole mass. The U T (p 2 ) are the transverse parts of the DR renormalized 
gauge-boson self-energies. Their explicit forms are given in Ref. [H], [17|. The form factors iT^ 



for one-loop corrected xtx.j G 11 vertices (G = 7, Z) have the following forms 
T^ = -r{v G LlJ P L + v%P R ) 

+F$ I ffP L + f$ r1 "p r + F° L (p 3 - p,yp L + F° R (p 3 - p,yp R 

~(vg f 5Zf fj + 5Z^vg l3 )YP L - \(vg f 5Z* + bZ^v G m ^P R . (6) 

The first line of Eq. @ contains the tree level couplings (|A.7|) with (e, gz-, U, V) in the DR scheme. 
Fy and F<f are the one-particle-irreducible (1PI) corrections to the vertices. Their explicit forms 
are given in the Appendix B. The last line gives the chargino wave function renormalization. 
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The wave function corrections 5Z L,R are determined in terms of the two-point function iKij(p) 
of charginos xl(~ P)xj(p) i n the DR mass basis. is decomposed as 

KM = ^(p 2 )j>P L + £g(p 2 )^ + E°(p 2 )P L + Ef*(p 2 )P R , (7) 

and SZ are then fixed by imposing well-known on-shell renormalization conditions for fermions 



18| . The diagonal parts of 5Z are 
5Zh= -Ef i (m 2 ) + — 



mi 



5Z« 



—m 



m,- 



- rrii 



mi 



Sf(m 2 ) + Er(m 2 ; 



(8) 



Xi 



is used for convenience. The term 



Here E'(p 2 ) = dE(p 2 )/dp 2 . The abbreviation m^ 

proportional to ImX^ in SZ^ comes from our convention to use real SZ R . In this paper we treat 
only cases where £^(m 2 ) is real (no CP violation). The off-diagonal terms (i 7^ j) are 



rm 



m 2 EfJm 2 ) + minijEl 



(ni- : ) + rriiEfjim 2 ) + m^*(mT 



SZ 



mt 



rriimjEfjim 2 ) + m 2 £^(m 2 ) + m^m 2 ) + m^*(m 2 ) 



'3 ij 1 



J' 



In addition, the pole masses of charginos are given by 

1 



mj(pole) 



rrii 



-mi 



+ E 



S^(m 2 ) + E^(m 2 ; 



(9) 



(10) 



In the corrected amplitude Eqs. (§, |5|), the renormalization scale dependence of the DR tree 
level parameters and that of the loop functions exactly cancel to O(a). However, the cancellation 
is quite complicated as we will see in the next subsection. 

In the numerical calculation we take the pole masses of gauge bosons (Z, W), and the standard 
model MS electromagnetic coupling a SM (M 2 ) as inputs. The DR gauge couplings are obtained from 
these parameters as discussed in Ref. fl3| . The chargino sector is fixed by giving pole masses of 
two charginos and tan/3(Mz). 

Finally, the spin-averaged differential cross section is written in terms of the amplitude as 



da 



dcos9 2s 167r 



h El-^l 2 - 



in 



Here YL denotes average over the initial electron and positron helicities and sum over the final 
chargino helicities. We use the helicity amplitude method [JTS] in the numerical calculation of the 
cross section. The relevant formulae are given in Appendix C. The phase space factor is given 
by 

1 



2(m._ + mz + )s + (m 2 



mt 



(12) 



s v Xj Xi 

Since a highly polarized electron beam will be available at future e + e~ linear colliders, in this 
paper we often present the cross section for an initial electron in a helicity eigenstate. Note that 
the chargino masses in /3% and in the wave functions (u(p 3 ), v(pi}) are the pole masses. 
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2.2 On-shell Renormalization of charginos 

The wave function renormalizations of charginos appear in the corrected amplitude Eqs. (01, 

|). They contain ultraviolet divergences from (E L , E R , T, D ) and, after DR renormalization, depend 
on the renormalization scale Q. This Q dependence cancels the implicit Q dependence of the DR 
mixing matrices (U, V), the gauge coupling g 2 in Eq. (||), and the explicit dependence of the 1PI 
vertex corrections in Eq. (ffi). However, this cancellation is quite complicated. For example, in 
Eq. (H) the Q dependence of the off-diagonal parts of 5Z L cancels both that of U in and that of 
Fy L . Moreover, the off-diagonal parts of 5Z L( - R - ) in Eq. @ superficially diverge when two chargino 
masses become degenerate. Therefore the forms of 5Z L ^ in Eqs. (^, |5|) can be inconvenient in 
real calculations. 

In this subsection, we reorganize the contribution of 5Z L ^ into a very convenient form, by 
utilizing the Q-independent effective chargino mixing matrices (U p , V p ). The loop contributions 
which compensate the running of (U, V) are then completely split from other corrections. 

We first notice that both the diagonal and off-diagonal parts of the chargino wave function 
renormalization 5Z L ^ can be implemented by making the following replacements in the couplings 
of the tree level amplitude, 

U}k —> U\ k + -UjjSZj^ , 

Vl - Vl + \v?5Zf k . (13) 

The corrections Eq. ([13|) are universal in any processes involving on-shell charginos. Remember 
that in Eq. (p~3| ) the mixing matrices (U, V) diagonalize the DR tree level mass matrix Eq. ( |A.2|) . 
The factors in Eq. (|T3"D come from the relations between the DR fields in the gauge eigenbasis 
the DR fields in the tree level mass eigenbasis xl , an d the on-shell renormalized fields x[ F i 

= UJjXJl = Ul^fffxkl , 
*jR = V SxJn = VS( zR f/kXkR ■ (14) 

We then introduce the effective mixing matrices of charginos, which are renormalization scale 
independent, and rewrite Eq. fll3| ) by these matrices. We first define the effective mass matrix 
M c {p 2 ) in the DR gauge basis as 

M c (p 2 ) = M C - t D {p 2 ) - l -M c t L {p 2 ) - \t R {p 2 )M c . (15) 

Eij are chargino two-point functions ipfipj in the gauge basis. They are related to £ in the DR 
mass basis as 

s L = ut L u ] , 

Y, R = V*t R V T , 

Mq is diagonalized by two unitary matrices, the effective mixing matrices U(p 2 ) and V(p 2 ), as 

M D {p 2 ) = V\p 2 )Mc{p 2 )u\p 2 ) , (17) 
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where Mr>{p 2 ) = diag(mj(p 2 )) is a real diagonal matrix. Note that (Mc, U, V, Md) are indepen- 
dent of the DR renormalization scale Q. 

We then give the forms of (U, V) and in terms of two-point functions S(p) of charginos. 
(U, V) are expanded as 

U{p 2 ) = U + 5U{p 2 ) = (1 + 5u{p 2 ))U , 

V(p 2 ) = V + 5V(p 2 ) = (1 + 5v(p 2 ))V . (18) 

Here 5u = 5U ■ and 5v = 5V ■ must be anti-hermitian from the unitarity of (77, V). Their 
diagonal elements must be then pure imaginary. 
The 0(a) expansion of Eq. flTT| ) gives 

M D (p 2 ) = (l + 5v*(p 2 ))V*M c (p 2 )U\l + 5u(p 2 )) 



M D - Z D (p 2 ) - l -M D H L {p 2 ) - h: R {p 2 )M D 
+5v*(p 2 )M D + M D 5u(p 2 ) . 



(19) 



The real parts of the diagonal elements of Eq. ( |T9| ) give m« = mj(pole) at p 2 = m 2 . The off-diagonal 
elements of Eq. (O) give the following relations 



5u i:j (p 2 



mf — wA 



Urn 2 + m 2 )^ + rmmjEZ + r^Eg + m^f ) (p 



R 



Hb 2 ) = " la foK 2 + m i) E 3 + "WjEy + + ™^*) (p 2 

for i j. Finally, the imaginary parts of diagonal elements give the relation 

1 



{-5u + 5v*)u(p z 



2m. 



D sr^D* 



By convention we set 5vu(p 2 ) = 0. The relations Eq. ( |Hp are then rewritten in terms of 



Uf^U^m 2 ), V^V^m 2 ), 



(20) 



(21) 



(22) 



as 



^iR 



N, 



1/2 



{U%{5 jk -\^ k {ml))Nl' 2 x k E 



rat 



R' 

kk ^kk 



m k 



^kk + ^kk ( m k) ■ 



D * 



(23) 



1/2 

N k is the real diagonal finite factor. The chargino wave function renormalization is then included 
by replacing U\ k and V^. in the couplings of the tree level amplitude by corresponding factors in 
Eq. pl), as 



4 - (U%(5 jk - l -^ k (rn 2 k ))N l k 12 



(24) 



The use of effective mixing matrices (U P ,V P ) has several nice features. First, the superficial 
singularity of SZ^j for degenerate masses is completely absorbed into (U p , V p ). We can then see 
that the original singularity just reflects the arbitrariness in the diagonalization of a matrix with 
degenerate eigenvalues. The absence of this singularity is similar to the procedure proposed in 
Ref. Second, the renormalization scale dependence in Eq. ( |25| ) only appears in H L ' R . The 



Q-dependent parts of the first equation of Eq. (j23|) become, up to O(a) 



1 ~ 

2 



(Sij - ^\ div )(U%Xk£ ■ (25) 



Eq. (p5|) takes the same form as the SU(2)xU(l) symmetric renormalization of ^>7". This property 
is very convenient both in theoretical considerations of the renormalization and in numerical 
calculations. Note, however, that {U p , V p ) are non-unitary at 0(a), unlike (U, V) and (U, V)(p 2 )f\. 

The wave function corrections in our process are expressed as follows. By applying the rule in 
Eq. (|24j), the s-channel form factors iT^ in Eq. (^]) are rewritten as 

rg M = -i»n] /2 n] I2 {v g Li3 p l + v%p r ) 

+ (1PI vertex corrections) , (26) 

where v^ R ) are obtained from i>£(m in Eq. ( |A.7| ) by replacing (U, V) by (U P ,V P ). The Q- 
dependence of the second line of Eq. ( PB"D exactly cancels that of the third line. The Q-independence 
of the rewritten form factor Eq. (pq) is thus more transparent than the original form Eq. (Jg). Sim- 
ilarly, the last factor in the t-channel amplitude of Eq. (£|) is rewritten as 



= V p *V p Nl' 2 Nf - liVtV^Um)) + Z$(m*)V p *V p ) . (27) 

In leaving, we comment that the effective matrix method given here can be applied to any 
process involving on-shell charginos, since the corrections of Eq. (^) are universal. This method 
can also be extended to other particles with flavor mixing, such as the neutralinos. 



2.3 Large limit 

We are interested in the limit where the squark mass is much larger than the masses of 
the charginos, the sneutrino, and the beam energy. Some corrections to the chargino production 
amplitude do not decouple in this limit but increase as logM^. This reflects the supersymmetry 
breaking in the effective field theory below the squark mass scale J|. 

*The effective matrices are unitary up to corrections of 0((m^2 — m^)/m^). 
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First, the effective chargino mass matrix Mc(p 2 ) receives non-decoupling corrections. For 
M~ 3> p 2 the effective mass matrix becomes independent of p 2 . The asymptotic form is obtained 
by replacing the elements of the tree level mass matrix Mc of Eq. ( A.2 ) by 



M, 



M 2 (Q) 



9g : 



1 + 



16tt 2 v 

Ky 2 + vl 



V2M W cos /3 
\/2M w sin 3 



167T 2 

V2M W cos 3(Q) 
V2M W sin P(Q) 




( -ir 



"7 =^ 



off 



, 5M W 



, 6M W 

1 + + 



M 



IV" 



cos /? 
5 sin /?' 
sin /3 



where 



167T 



16tt 



16tt 



, 5 cos /3 
cos /3 

2 (5 sin /3 
sin /3 

,5M W 



sin 2 /3 



M 2 Ml 
% ln "T^T - ln 



" 2/3 

•-cos p 
2 P 



Mi 




Q 2 



3 M l 
yi cos 2 8 ln — | 



Vb In 



14' 

o M~ 
2 y * M M 2 , 



9 2 



9 2 
+ 4^ 



(28) 
(29) 
(30) 
(31) 

(32) 
(33) 
(34) 



(3R - 2) \n(R - 1) + irln 



R-l 
R 



do not decouple in the large Mq limit. This effect is very 



with R = m 2 /M^. The corrections to the diagonal elements can be absorbed into the effective 
mass parameters M| ff and /z eff and are not interesting within the context of the MSSM. By 
contrast, the corrections to the gaugino-Higgsino mixing masses cannot be absorbed into unknown 
parameters such as tan/3. The squark loop corrections to the effective mass matrix Mc, and 
effective mixing matrices (U p , V p 

important if the gaugino-Higgsino mixing is not highly suppressed 

In the s-channel amplitude all other squark loop corrections decouple in this limit. The squark 
loop corrections from the gauge boson self energies I1 T decouple after the gauge couplings are 
renormalized. The factor N l l 2 in Eq. ([23]) approaches to 1 in this limit. Finally, the non-decoupling 
terms in E L ^^ in Eq. (|26| ) exactly cancel the Fy L ^ terms of the 1PI vertex corrections in Eq. (j^j. 
This result is consistent with the universality of gauge boson interactions. 

By contrast, the 0(logMg) terms in T, L ( RS) remain in the t-channel amplitude and are very 
important. This is the origin of the "super-oblique corrections" discussed in Refs. [[11], [T2|, [T3"|j . 
We point out that the corrected t-channel amplitude takes a very simple form for sufficiently 
heavy squarks. In this case, the corrected amplitude is obtained from the tree level one by the 
replacement 



gi(Q)v^v 3 



9, 



evW 



(35) 
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Here we used the fact that Sy,(m|) [f ^ 1) is insignificant for sufficiently heavy squarks. The 
parameter g e ^, which is renormalization scale independent, is interpreted as the effective evW 
coupling. g e ~yr deviates from the corresponding gauge coupling, gf M (Q). Its asymptotic form is 



Here we note that the dependence of g e -^(m^) decouples in the large Mg limit. 

3 Numerical results 

In this section we describe the dependence of the chargino production cross section cr(e~e + — > 
X1X1) on various MSSM parameters. The production cross section is a function of the gauge 
couplings, the Z-boson and sneutrino masses, and the chargino masses and mixing matrices (pa- 
rameterized by M 2 , fi, M\y and tan/?). In the proposed colliders the electron beam can be highly 
polarized, therefore we often show the production cross section with a polarized electron beam. 
We denote the cross section as (Tl{r) when the initial state electron is left handed (right handed). 

In the gaugino region (M 2 <C Xi is wino-like, and the amplitude receives both t-channel 
and s-channel contributions, unless the initial electron is right handed. If the electron is right 
handed, the t-channel amplitude vanishes because of the absence of a Wcrv coupling. In the 
opposite limit, <C M 2 , the lightest chargino is Higgsino-like. Since the Higgsino couplings to 
the first and second generation (s)leptons are negligible, in the Higgsino limit only the s-channel 
amplitude contributes. Finally, when M 2 ~ both charginos have large gaugino and Higgsino 
components, and they are somewhat degenerate in mass. In this region of parameter space the 
chargino mixing matrices relevant in the production cross section are sensitive functions of tan/?, 
which enters in the off-diagonal elements of the chargino mass matrix. 

Formulas for the one-loop corrected chargino production cross section are given in the previous 
section and the Appendix B, including quark and squark loop effects. The t-channel amplitude 
depends on the effective coupling g.~^, the effective chargino mixing matrices V p , and decoupling 
corrections. The s-channel amplitude depends on the usual gauge couplings, the effective mix- 
ing matrices U p , V p , and decoupling corrections such as the 1PI gauge-chargino-chargino vertex 
correction. 

Both amplitudes depend on squark masses m^, squark mixing angles 9^, and quark-squark- 
gaugino(Higgsino) couplings. In this section, we present our results assuming a universal soft 
breaking squark mass M-x and a universal trilinear coupling A at the weak scale. These parameters, 
along with \i and tan/3, determine the squark masses and mixing angles. The third generation 
quark- squark- Higgsino couplings depend on the top and bottom Yukawa couplings y t and y^. As 
shown in Eqs. (|3Ti|-pM|) , the heavy top quark can give rise to a sizable correction proportional to 
yflnM-, which enters in the off-diagonal elements of the effective chargino mass matrix. The 
Yukawa couplings are also involved in the 1PI vertex corrections when the final state chargino is 
Higgsino-like. The top Yukawa coupling is very large when tan /3 — > 1 while y\> is substantial when 
tan/3 > 30. 



[0 0] 




(36) 
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In the following we will consider XiXi production in the three cases where the lightest chargino 
is predominantly gaugino, predominantly Higgsino, and in the mixed region. We will refer to the 
parameter sets listed in Table We now discuss the three regions in turn. 



name 


description 




m x~ 2 


tan (3{M Z ) 


vriv 


A 




sgn n 


Gl 


gaugino region, > M 2 


200 


800 


2 


100 





500 


-1 


HI 


Higgsino region, < M 2 


200 


800 


2 


100 





500 


-1 


G2 


gaugino region, > M 2 


172 


512 


4 


240 





500 


-1 


H2 


Higgsino region, < M 2 


172 


512 


4 


400 





500 


-1 


M 


mixed region, \fi\ < M 2 


172 


255 


4 


240 





500 


-1 



Table 1: Five parameter sets. All entries with mass units are in GeV. 



3.1 Gaugino region 

In Fig. |l](a) we plot the chargino production cross section a(e~ L e + — > XiXi) versus (solid line), 
for the Gl gaugino region parameter set of Table |l]. In the gaugino (or Higgsino) region the diagonal 
elements of the effective chargino mass matrix Mc are fixed by the input chargino masses, so they 
are independent of M~. Conversely, the DR parameters M 2 (M 2 ) and fi{ii) vary as Mq increases. 
The effective mixing matrices U p , V p contain non-decoupling log(M^) + constant corrections. 
These corrections arise from corrections to the off-diagonal elements of the effective chargino mass 
matrix given in Eqs. (|30H 34l). They contribute in both the s- and t-channel amplitudes. However, 
the dependence of the mixing matrices on the off-diagonal elements of the effective chargino mass 
matrix is suppressed for this set of parameters, so V£, ~ 1 over the whole range of squark 
mass shown. The positive correction proportional to log M-x in Fig. |I](a) is therefore primarily due 
to the loop correction to the effective coupling g e -^- 

The remaining corrections vanish in the large limit. These remaining corrections can be 
divided up into oblique and non-oblique parts, each of which satisfies decoupling. The non-oblique 
part consists of the 1PI vertex correction and the associated chargino wave function renormaliza- 
tion. In the following when we refer to the vertex correction we mean this combination. The 
vertex correction is somewhat complicated, so it is worthwhile checking whether this gauge and 
scale invariant correction can be neglected. The cross section calculated without including the 
vertex correction is shown by the dotted line in Fig. |l|(a), and the ratio between the cross section 
without the vertex correction and the full one-loop cross section, (j™~ vtx /o^, is shown in Fig. |l](b). 
The maximum effect of the vertex correction is less than 0.5% of the total cross section for this 
choice of parameters. The vertex correction is negligible compared to the sensitivity to in 
future experiments. 

As seen in Fig. |l](a), the left-handed cross section ol increases by about 14% as varies 
from 300 GeV to 3 TeV. The sensitivity to M-x depends on mc, and as discussed in Ref. Jll 
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Figure 1: (a) The one-loop chargino production cross section ai as the function of the soft squark mass 
M~ for the gaugino-like parameter set Gl of Table 1 (solid line). The positive correction proportional 
to logM~ is due to the loop correction to q ~~. The dotted line shows the cross section without the 
Z{l)x~X + vertex corrections, (b) The ratio between the cross section without the gauge vertex corrections 
and the full one-loop cross section for the same set of parameters. The vertex correction is less than 0.5% 
of the total cross section. 

° L (pb) 



4000 

> 
o 

S 1000 
400 

100 150 200 250 300 350 400 450 500 
nyGeV) 

Figure 2: Contours of constant ox in the (M~,mp) plane for the parameter set G2 of Table 1. ol 
increases (decreases) with M~ if mp 1z 200 GeV 300 GeV) . 




In Fig. fj we show contours of constant in the (mp, M~) plane, for the G2 parameter set of 
Table |l[ As increases, increases if mp is less than 200 GeV, while it decreases if mp is 
greater than 300 GeV. For mp ~ 250 GeV, gl becomes insensitive to M-x. The dependence on 
Mq from the t-channel amplitude is negligible in the limit mp ^> y/s since the t-channel amplitude 
scales as 



3.2 Higgsino region 

In Fig. 13 we show the Mp. dependence of <jl when the chargino is Higgsino-like. We take parameter 
set HI of Table |I|. The diagonal elements of the effective chargino mass matrix Mq are fixed by 
fixing the chargino masses. As in the gaugino region, the mixing is suppressed, | Vi_2 1 5 I C/12 1 — 1- The 
one-loop cross section including (not including) the vertex correction is shown by the solid (dotted) 
line. The cross section changes by less than 0.5% as varies from 300 GeV to 3 TeV. Such a 
weak Mq dependence in Higgsino-like chargino production is expected from our observations in 
Sec. 2. Although the large top quark Yukawa coupling is involved, the vertex correction remains 
small, less than 1%. 
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Figure 3: The cross section gl with/without the gauge vertex correction vs. M~ for the parameter set 
HI of Table 1 (solid/dotted). The lightest chargino is Higgsino-like. Both the dependence on M~, and 
the effect of the vertex correction is very weak. 
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Figure 4: <tl vs. M~ for the mixed region parameter set M of Table 1 (solid). The s-channel and the 
t-channel cross sections are shown by the long dashed and short dashed lines. The cross section receives 
important corrections from the logM-x dependence of U p and V p (see text). 



3.3 Mixed region 

In the mixed region (M2 ~ |/i|), the full one-loop cross section receives important corrections 
proportional to logMg through the corrections to the effective mixing matrices U p , V p , as well 
as logMg corrections from the effective coupling g e -y^- We illustrate this in Fig. ^, which shows 
the production cross section for parameter set M of Table [I]. For this choice of parameters M 2 
and \fi\ are both near 200 GeV, so the chargino mass eigenstates are fully mixed (|^n| 2 ~ 0.6). 
In the figure, <7l increases by 4% as varies from 1 to 10 TeV (solid line). The destructive 
interference between the t-channel and s-channel amplitudes accounts for this insensitivity. The 
22% reduction in the t-channel cross section (short-dashed line) is due to an 8% reduction of V£ 
and a 2% increase of g ej> $r- The s-channel cross section depends on both U p and V p , and decreases 
by 7% (long-dashed line). 



3.4 Comparison of Mg and tan/3 dependencies 

We now compare the and tan/3 dependence of the chargino production cross section. In 
Figs. ||(a)-(e), we show contours of constant cross section in the (Mq, tan/3) plane. In Fig. |](a) 
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Figure 5: Contours of <tl (Figs, (a), (b) and (d)) and ctr (Figs, (c) and (e)) in the (M~, tan/?) plane. The 
parameter sets from Table 1 are (a) G2 (gaugino-like) with yfs = 400 GeV; (b) and (c) H2 (Higgsino-like) ; 
(d) and (e) M (mixed) with mr, = 400 GeV. The M~ dependence is strong in Fig. (e), though the tan/3 
dependence is more prominent. The cross hatched regions are excluded by the chargino mass constraints. 



we use parameter set G2 of Table [I], with a/s = 400 GeVf]. Since Xi is gaugino-like, <jl depends 
on my. We see that <jl is insensitive to tan/3 in this case. It is almost constant in tan/3 when 
tan/3 > 5. On the other hand, decreases by 10% when changes from 300 GeV to 10 TeV 
due to the correction to g e -^- 

In Figs. |5](b) and (c), we show the results for Higgsino-like Xii using parameter set H2 of 
Table [I]. Over the entire variation of tan/3 and Mg considered, varies by less than 1.5% 
(Fig. §(b)). The cross section or is relatively more sensitive to tan/3 (Fig. 0(c)), but the absolute 
change in the cross section is smaller than in the <jl plot. In both cases the dependence is 
very weak because of the very small mixing, and the absence of a t-channel coupling. 

In the case of large gaugino-Higgsino mixing, the cross section is more sensitive to tan /3 than 
to Mg. In Figs. 0(d) and (e) we show the chargino production cross section in the mixed region, 
with parameter set M of Table 0, except rrty = 400 GeV. The cross hatched region tan (3 Z 4 
shown in these plots is excluded because it is not possible to obtain the specified chargino masses 
in this region. We find a strong dependence on tan/3 for both a L (Fig. 0(d)) and or (Fig. 0(e)). 
The dependence is very small in Fig. (d) due to the interference between the s- and t-channel 
amplitudes. In general it can be large. For example, in Fig. (e) at tan {3 = 4, or changes from 45 
fb to 62 fb as changes from 300 GeV to 10 TeV. 

The large dependence of the one-loop cross section on in the mixed region is caused by 
the strong sensitivity of the cross section to the off-diagonal elements of the effective chargino 
mass matrix. In Ref. J7j, it was claimed that a 4% measurement of cxr results in the constraint 
3.9 < tan/3 < 4.1 (at tan/3 = 4). The tree level chargino masses were fixed in their determination. 
We see from Fig. 0(e) that for a given value of <7r the central value of tan/3(M^) can shift by 0.5, 

tThe parameters are those used in Ref. [jfl]. We take -Js = 400 GeV for Fig. |^(a) because of the accidental 
insensitivity of <jl to M~ at ^/s = 500 GeV. 



14 



°k(fl>) 




Figure 6: Contours of ctr in the (A/M~, M~) plane for the mixed region parameter set M of Table 1. 
The A/M~ dependence is larger for smaller M~. 

depending on Mq. 

We find that while the vertex correction is not substantial compared to the experimental 
sensitivity to the cross section, the corrections due to U p and V p can be. For example, the 
dependence in Fig. |(e) is almost entirely due to U p and V p . At tan/3 = 4, <tr changes by about 
37% as Mq varies from 300 GeV to 10 TeV. In contrast, the vertex correction is less than 1.4% 
over this range. 

3.5 Squark mixing effects 

Left-right squark mixing effects give rise to important corrections in the mixed region. The stop 
and sbottom mixing angles are controlled by A, fi, and tan j3 as described in Appendix A. In Fig. [] 
we show contours of constant ctr in the (A/Mp:, M~) plane for the chargino in the mixed region. 
We use parameter set M of Table [y. The cross hatched region of Fig. |] is excluded either because 
of the chargino mass constraint or because m~ < 0. The cross section shows strong dependence on 
A/Mq when is small. For example, when A/M^ varies from —2 to +2 with = 345 GeV, 
the cross section changes from 61 fb to 39 fb. The stop mass eigenstates become fully mixed at 
large \A\/Mq. The top squarks are well split when A/M^ = -2 (m~ tl = 128 GeV, m~ t2 = 528 GeV), 
while at a point of small mixing, A/M^ = 0.12, the top squarks are nearly degenerate (m^ = 383 
GeV, mi 2 = 386 GeV). The A dependence of a R decouples at large as 1/Mg for fixed A/Mq. 

If only the two chargino masses and <jr are measured it may not be possible to disentangle 
the dependence of the cross section on squark mixing from the dependence on tan /3 and M^. For 
example, at = 345 GeV in Fig. |6], we see the cross section is 61 fb at A/M^ = —2. We can 
find the same cross section with the same chargino masses at A = by changing tan (3 from 4 to 
3.6. It may be necessary to measure A from other quantities. The stop masses and mixing angle 
can be constrained if 2 and a(e~e + — > tit\) are measured JEJ. Combining these measurements 
with measurements of [i and tan j3 from other processes, A t can then be determined. 

Because left-right squark mixing arises from SU(2) x £7(1) gauge symmetry breaking, it con- 
tributes to the violations of the relations between the tree level chargino and neutralino masses. 
We can utilize this dependence in efforts to constrain the values of A and tan/3. In Fig. [?] we 
show contours of A = varying tan/3 and of tan/3 = 4 varying A in the (wi^o, <7r) plane (solid 
lines) . We fix Mq = 345 GeV, and the other input parameters are fixed as in Fig. ||. The contours 
terminating at m^o = 201.4 GeV have tan/3 = 4 while the contours terminating at m^o = 212 A 



15 



0.12 

0.10 

0.08 

Igj 0.06 
D 

0.04 
0.02 
0.00 



A or tan(5 dependence in mixed region 



190 



170 




200 210 

m„„(GeV) 



160 



Figure 7: Contours of A = (tan/3 = 4) varying tan/3 (A) in the (m^o, <tr) plane (solid lines), with 
M~ = 345 GeV. The variation of m^o with m^o for the same set of parameters is also shown (dashed 
lines). 



GeV have A = 0. For values of tan/3 slightly above 4, we cannot find solutions with the given 
chargino masses. For a given <jr, we find m^o differs up to 3 GeV between two contours. We also 
show the variation of m^o with m^o when A or tan/3 is varied (dashed lines). For fixed m^o, m^o 
differs up to 2 GeV between the two curves. If we can measure the chargino and neutralino masses 
within 1 GeV or less, it could help to single out the effect of squark mixing |. 

Notice the squark mixing dependence of the radiative correction is mainly due to the correction 
to U p and V p , not the vertex correction. For example, at M-x = 350 GeV, <jr changes from 47.3 fb 
to 41.0 fb as A/Mq changes from to 2, while the cross section without the vertex correction 
varies from 46.7 fb to 40.0 fb, a 1.4% to 2.4% effect. Although it is larger than the < 1% effect 
found in the gaugino and Higgsino dominant regions, it is unimportant compared to the strong 
dependence of the cross section on A, tan (3, and Mx. 

3.6 Comparison and approximations 



We should briefly comment on the comparison of our results with those of Ref. [15"]. The results in 
Ref. |HJ are obtained by including top, stop, bottom and sbottom loops only. They underestimate 
the ^|log(M^) corrections to g e -^/g2 and U p , V p , which depend equally on the 1st and 2nd 
generation (s) quarks as the third generation. We have checked the t-channel exchange of the 
sneutrino has a substantial effect on the total cross section in some of their plots, and including 
the contributions of the first two generations significantly alters the results. 

Our comparisons with their results show large numerical differences. For example, for the 
parameters corresponding to their Fig. 4 at tan (3 = 0.5, we find the one-loop cross section decreases 
by 1.2% as Mg = Af\ varies from 200 to 1000 GeV. Their results show a 17% increase in the cross 
section. Notice they take the gaugino mass parameter M 2 (Mz) as input. Taking this unphysical 

•fAn excellent measurement of the chargino and neutralino masses may be achieved at proposed colliders. 
A recent study shows that it should be possible to measure the lighter chargino mass with an accuracy of 30 to 300 
MeV by measuring the cross section in the threshold region [^3| . 

*We refer to their A. We use the opposite sign convention for A. 
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mass parameter as input generally leads to larger M~ dependence. In the example just mentioned, 
taking m~- as input reduces the M-x dependence from 1.2% to 0.4%. 

A 2 V 

We also find smaller differences between the tree level and one-loop cross sections. However, 
this is not surprising, since the definition of the tree level cross section is somewhat arbitrary. Our 
tree level cross section is determined by the two chargino masses, M w , M z , m„, tan/3, and the 
effective theory (i.e. standard model) MS gauge couplings. The tree level cross section depends 
on the choice of the scale of the effective theory gauge couplings. An appropriate scale is found 
by considering the Higgsino production cross section in the limit Mz -C y/s <C M^. In this 
limit the squarks are completely decoupled. Because of the constant correction in the quark 
loop oblique correction, the cross section in the effective theory is equal to the cross section in 
the full theory if the effective theory gauge couplings are evaluated at the renormalization scale 
Q = exp(— 5/6) a/s ~ y/s/2. Hence, our tree level cross section is evaluated with effective theory 
gauge couplings evaluated at the scale y/s/2. With this choice, the tree level and full one-loop 
cross sections are nearly equal when the squark corrections decouple. 

We have already discussed that, for practical purposes, it is safe to neglect the vertex correction. 
We will now consider two approximations to the remaining corrections. In the first, the effective 
theory approximation (ETA), we use the effective coupling g e -^ and the effective mixing matrices, 
U p , V p . In this approximation some 1/M~ corrections are included. In the second approximation, 
the "log + constant" approximation (LCA), we strictly keep only the non-decoupling squark 
corrections, i.e. we include only the corrections of the form log(M^) + constant. The effective 
coupling in the log + constant approximation, (Jlca, is given in Eq. (0). The effective mixing 
matrices in the LCA are found as follows. The off-diagonal elements of the effective chargino 
mass matrix in the LCA are given in Eqs. (|30}-|34|). The effective mixing matrices in the log 
+ constant approximation, U£ CA , V^ CA , are then determined from the two chargino masses and 
these off-diagonal elements. Notice that <?lca, U£ ca , V^ ca , and the effective theory couplings, are 
renormalization scale independent. 

In Fig. |3| we show the ratio of the various approximations to the full cross section, versus M^. 
We plot the ratio of unpolarized cross sections in the gaugino, Higgsino, and mixed regions in 
Figs. |K a), (b), and (c), with parameter sets Gl, HI and M of Table [l], respectively. The cross 
section without the vertex correction is shown by the dotted line. The ETA result is shown with 
the dot-dashed line, and the LCA result is indicated by the dashed line. 

There are two factors which contribute to the deviations from unity in the large region 
in the ETA and LCA results. For one, these approximations are calculated with effective theory 
couplings, while the full calculation is calculated with full theory couplings. This mismatch causes 
discrepancies of order (alogM-) 2 . These discrepancies give some indication of the expected 
magnitude of two-loop corrections. Another reason why the approximations can disagree in the 
large Mq limit is that the scale used to evaluate the s-channel tree level gauge couplings is y/s/2. 
The scale which should be used to get exact agreement in the decoupled regime is somewhat 
different, depending on M z , m t , and y/s. 

In all three figures the vertex correction is less than 1%, so the "no vertex" approximation is 
a good one, even for very small squark masses. The ETA also works well, better than 1% except 
at M-x < y/s/2 in Fig. (b). The LCA works as well as the ETA, except in the gaugino region with 
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Figure 8: Various approximations to the cross section divided by the full one-loop cross section. The 
results for the "no vertex" approximation, effective theory approximation (ETA), and "log + constant" 
approximation (LCA) are shown. Figures (a), (b) and (c) show results in the gaugino, Higgsino, and 
mixed regions of parameter space, respectively. (Parameter sets Gl, HI, and M of Table 1 are used.) 



Mg < a/s and in the mixed region with < 1.5-y/s. By comparing the ETA and LCA results we 



Q 



see that the 1/M~ corrections included in g e -y^, U p and V p can be essential in obtaining a good 
approximation, even for squark masses as large as 1.5y/s. 



4 Uncertainty in the chargino production cross section 
measurement 

4.1 Previous analyses 

In this section we revisit previous studies of chargino production and decay [|7|, [Tl[. Chargino 
production can be studied in e + e~ collider experiments by observing their decay into v£\\ or 
qq'x\i with signals £2j + missing momentum or 4j + missing momentum. 

In Refs. [7L |T]]] the probe of the supersymmetric relation of SU{2) gauge/gaugino couplings 
g 2 M =g ej> w was considered based on the MC study of the £2j mode at the point in parameter space 
{fi, M 2 , tan (3, M 1 /M 2 ,m c ,) = (-500 GeV, 170 GeV, 4, 0.5, 400 GeV), where m-+ = 172 GeV and 
86 GeV. The analysis of Ref. |7j assumes y/s = 500 GeV and that no direct production of v 



1 

is available. This results in a poor constraint on g ej> ^- In Ref. |11| the authors consider the same 
point in parameter space, except they assume is measured directly, and yfs can be tuned. 
Further, they assume that the uncertainty of the theoretical input parameters in the chargino 
and neutralino sector, the acceptance of £2j events, and the dependence of the acceptance on the 
theoretical input parameters are independent of and y/s. Under these assumptions, they find 
Sg ~wl gt M = 2% with m„ = 240 GeV and a/s = 400 GeV. This result is considerably poorer than 



the result 8g e -^/ g 2 < 0.6%, estimated in the sneutrino production study of Ref. |T3| . 

The purpose of this section is to provide a critical discussion of the analysis of Refs. [0, [TlH . 
We point out that the poor constraint found in Ref. 0] results from the low acceptance of the 
£2j mode (found in Ref. 0), and the low acceptance is further traced back to the (special) choice 
of parameters. We also point out that the signal acceptance and dependence of the acceptance on 
the theoretical input parameters can be strongly dependent on a/s. Therefore the estimate given 
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in Ref. []TT] cannot be trusted at other values of y/s unless a dedicated MC simulation is provided 
for both the signal and background. The bottom line is that the constraint on g e ~^ can be greatly 
improved for more generic parameters, and by optimizing the beam energy and cuts. 

Before going into the details of their simulation, we shall discuss the background to the £2j 
signal for the case where xt decays exclusively into y^W as follows 

e + e~ -> xtXi 

- W+W'xlxi 

-> Ivqq' + i> T ■ 

We assume Xi is the stable LSP, so it escapes detection and (along with the neutrino) gives rise 
to missing momentum in these events. 

This process suffers from W^-boson pair production background. In the background events the 
total momentum of the W-boson pair is balanced in the transverse direction, but the observed 
transverse momentum is not balanced, due to the escaping neutrino. Hence, the discrimination 
between the signal and background is difficult. 

In the MC study of Ref. @, the following cuts are made to reduce the background from 
PF-boson pair production in the £2j mode: 

a) existence of an isolated hard lepton; Eg > 5 GeV, 9 g i > 60 ° 

b) i> T > 35 GeV 

C) #acop >30° 

d) m £l/ISR > 120 GeV 

e) -Qecos6 had , Qtcos9 e < 0.707 

Cuts b), c) and d) are set to reduce the W-pair events produced nearly back to back in the 
transverse direction, while keeping the supersymmetric signal. The cut e) is designed to remove 
the large forward peak of the WW events. 

Although these cuts are standard ones to improve the signal to noise ratio, the acceptance of 
the signal turns out to be small, 

^ = —^- = 11.9%, (37) 
a L BiB h L 

resulting in S/N = 1 at the previously mentioned point in parameter space. 

Our knowledge of the acceptance is limited by the errors of the underlying parameters (//, 
Mi, M 2 , tan/3, mp). The systematic errors on 77 are estimated as A^ sys = 0.55% by allowing 
the underlying parameters to vary so that itl^+, m ^+? an d m^o vary within 2 GeV of their input 
values. 

While A^ sys itself is small, the error in the cross section due to the acceptance uncertainty 
turns out to be large, i.e. Arj/r] = 5%. This is comparable to the change in the cross section when 
nig changes from 1 TeV to 4.5 TeV with = 240 GeV and yfs = 400 GeV. 
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No question was raised concerning the small acceptance of the £2j mode in Refs. ]7|, [IT]], but 
the value should be contrasted with other MC studies, which generally claim acceptances of 30 to 
50% for SUSY signals [|, ^ |2H| . The low acceptance in Eq. flUTFD is a consequence of the special 
choice of parameters. At the point in question, m~+ = 172 GeV, and m$o = 86 GeV, so the mass 
difference between the parent and the daughter particles Am = m~+ — m^o — mw is only 6 GeV. 
In the rest frame of xf, both the PU-boson and Xi are nonrelativistic. When the parent xf is 
boosted, the angle between the xf momentum and the VF-boson, and also the momentum spread 
of the W, can be very small. When the charginos are produced at y/s = 500 GeV, < 15°, 

and the VU-boson pair momenta and Xi momenta are roughly balanced. 

Notice in the cuts listed above we are relying on large O^+w an d large missing transverse 
momentum to separate the signal from the WW background, but neither of these attributes are 
characteristic of the signal events. The small acceptance merely results from the fact that the W 
pair from the signal and the background have very similar kinematics at this particular point in 
parameter space. 

The small acceptance has a direct effect on the acceptance uncertainty. The signal event 
distribution in the ($t, acop , m&,) space sits near the background distribution and therefore near 
the cut region. When the input parameters are changed slightly within their error, the signal 
region also changes in the {^T,&a,co P , m iu) space. Because the accepted number of events for the 
input parameters is so small compared to the total number of reconstructed IV-pair events, a 
small change in the signal region easily changes the acceptance by several percent. Relatedly, 
0^+ w is a rather sensitive function of Am when Am is small. Note that the systematic error of 
the acceptance is estimated by changing Am from 2 to 10 GeV in Ref. 0. 

To illustrate the kinematics, we show the acoplanarity angle distribution of W pairs recon- 
structed from e + e~ — > xfx.i~ y WW XiXi ~~ > 4j + ]/>t events in Fig. Q. The 4j mode also suffers 
from the WW background, however the cuts to remove the SM background are far simpler than 
those of the £2j mode[[ 

Fig. |] shows the acoplanarity angle distribution after applying the cuts to reject background 
W-pairs given in Ref. ||, except for the acoplanarity angle cut # acop > 30°. To generate MC 
events, we modified the event generator of Ref. M, and we used the JLC detector simulator ||. 
The effect of initial state radiation is included. The distribution shown by points with error bars 
corresponds to our standard input parameters (/i, Mi, M 2 , tan/3, mp) = (—500 GeV, 84.6 GeV, 
170 GeV, 2, 400 GeV), resulting in m~+ = 176.6 GeV and m^o = 86.9 GeV. The parameters are 
chosen so that Am = 9.4 GeV, larger than that in Ref. @. The size of the error bars and the 
central values are determined from 10000 generated events, corresponding to / Cdt = 16 fb _1 . On 
this same plot we also show two distributions corresponding to Mi = 90.6 GeV (short-dashed) and 
Mi = 78.6 GeV (long-dashed). These distributions correspond to m^o = 92.2 GeV and m^o = 80.3 
GeV, respectively. The difference in Am for these two curves is 12 GeV, large enough to create 
a statistically significant difference in the event distribution, given the somewhat small integrated 
luminosity we are considering. These curves are normalized to have equal numbers of eventsf]. 

^The 4j mode suffers from the SUSY background due to e + e~ — > X2X2 ~^ ^JXiXi, which may be hard to 
distinguish from the chargino signal. We discuss the 4j mode here for more or less illustrative purposes, although 
it is possible to extract more physics information by including this mode in a combined fit. 

*The total number of reconstructed WW events depends on the chargino/neutralino mass differences. For 
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Figure 9: The number of accepted chargino production events vs. the acoplanarity angle, with yfs = 500 
GeV and / Cdt = 16 fb _1 . The points with error bars are for Mi = 84.6 GeV, and the long-dashed 
(short-dashed) line corresponds to M\ = 78.6 (90.6) GeV. See text for other parameters. 



The acceptance changes drastically for the different cases. Implementing 6 acop > 30°, we find 
the acceptance varies by a factor of two. The acceptance is correlated with Am, which determines 
the maximal x + W angle. Below we list Am, the maximal angle 0™+^, and the acceptance for the 
three cases. 

m^o Am 8^?^ rj 

92.2 GeV 4.1 GeV 13.7° 14.1% 
86.9 GeV 9.4 GeV 20.7° 19.3% 

80.3 GeV 16.0 GeV 27.2° 27.3% 

Varying m^o by 2%, we expect about a 2% change in the acceptance. This corresponds to A?]/?] 
of about 10%. 

Notice we require 8 acop > 30° to reduce the SM background, while the maximal 8 acop is 28™+^ 
if the reconstructed jet momenta are identified with the quark momenta. For the sample with 
Am = 4.1 GeV, the events are accepted by virtue of the finite resolution of the jet axis. As 8^+w 
increases above half of the acoplanarity angle cut, the accepted number of events increases linearly 
with 8^+ w . On the other hand, if Am is so large that #~+^ 3> 8 acop , then most of the events pass 
the cut by a wide margin. In particular, most events are accepted regardless of several percent 
variations in the input parameters. Therefore, the acceptance error is much smaller in a generic 
region of parameter space. We expect the acceptance uncertainty to scale roughly inversely with 
the acceptance, for sufficiently large Am. We will examine this conjecture later by an explicit 
example. The uncertainty itself depends on the mode under consideration and the cuts applied, 
as we discuss in the next subsection. 

For the £2j mode the situation is less clear. Each of the three cuts, the $t cut, the acoplanarity 
angle cut, and the m&, cut, causes roughly the same reduction in the number of signal events. 

example, rejection of the forward going jets (W's) gives such dependences. However, as discussed below, this mass 
sensitivity is small compared to the uncertainty in the acceptance due to the acoplanarity angle cut. 
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Because of the missing momentum from the escaping neutrino, each cut yields smaller reductions 
compared to the Aj mode. However, each acceptance dominantly depends on the parameter Am. 
In particular, the acceptance is larger for larger Am for all of these cuts. When Am is small the 
signal region significantly overlaps the background region. In that with the Aj mode, the 

variation of Am within the error is largely responsible for causing events to move into and out of 
the accepted region. 



4.2 Improving the measurement 



We found in the previous subsection that the point in parameter space considered in Refs. [0, [TT 
must be regarded as a pessimistic case. At other points the acceptance will generically increase, 
leading to a decrease in the acceptance uncertainty. However, nature may equally well choose any 
point, so let us reconsider this point for a moment, and seek a procedure to reduce the acceptance 
uncertainty. 

One possibility is to reduce the chargino and neutralino mass errors, especially the error on 
their mass difference. When there is no correlation between Sm-+ and Sm^o, the largest (smallest) 
acceptance comes from the point where m-+ — m^o becomes maximum (minimum) within the 
mass errors. For example, in the Aj mode, when we change M\ and M 2 so that (Am~+, Am^o) = 
(+2 GeV,— 2 GeV), we find the acceptance increases by 6%. With (Am^+, Am^o) =(+2 GeV,+2 
GeV) we find the acceptance increases by 1.7%|. Hence, the acceptance is over three times more 
sensitive to absolute changes in the mass difference than the mass sum. 

In the standard technique to determine a particle's mass from the energy distribution of one of 
the daughter particles in two body decay, the mass difference between the parent and the daughter 
particle is measured better than the individual masses, especially when the parent particle is 
significantly boosted. (See examples in Refs. P, p4fl.) So, at generic points the uncertainty in the 
acceptance is much smaller than one would expect from uncorrelated mass errors. However, in this 
particular example, the acceptance is smaller near the endpoints of the energy distribution. This 
is because the daughter particle has a maximal energy when it goes in the same direction as the 
parent particle. In such a case, the acoplanarity of the event comes only from the other chargino, 
leading to small statistics near the endpoints. Therefore, we expect that the energy distribution 
is less sensitive to the chargino/neutralino mass difference for the case given in Ref. |7|]. 

The uncertainty of the acceptance may be reduced by increasing the acceptance itself. One can 
increase the acceptance easily in the Aj mode by reducing yfs. In Fig. |H] we see the acoplanarity 
angle distribution of ly-pairs is much flatter for a/s = 400 GeV. The angle #™+^ = 41.3°, so a 
large number of events pass the cut > 30°. We find the acceptance increases to 54.8% from 
19.3% at our standard point (where m^+ = 176.6 GeV and m^o = 86.9 GeV). The acceptance 
increases by 2.4% with m^+ = 176.6 + 2 GeV and m^o = 86.9 — 2 GeV. The uncertainty in the 
acceptance from the chargino/neutralino mass errors is therefore reduced by factor of 7 relative 
to the v/i = 500 GeV case0. 

§At (+6 GeV, +6 GeV) we find the acceptance is 24.3%. We then assume a linear dependence on the Am^ to 
obtain the estimate for (+2 GeV, +2 GeV). 

^However, we find the acceptance error including the WW reconstruction efficiency is only reduced by a factor of 
2.6. The dependence of the VF-pair reconstruction efficiency on the chargino/neutralino masses comes in through 
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Acoplanarity angle distribution of reconstructed W pairs 




Figure 10: The acoplanarity angle distribution of M^-pair events from chargino production, at = 400 
GeV. See the text for other parameters and cuts. 



For the £2j mode, reducing y/s may not result in a larger acceptance. For this mode, many cuts 
are needed to reduce backgrounds and the relevant distributions have different ^/s dependencies. 
For example, the $t > 35 GeV cut rejects more events at smaller y/i, because the narrower 
allowed range of Xi momentum leads to observed events which are more balanced in transverse 
momentum. Since the WW background also has a softer ]/>? at smaller y/s, it may be beneficial to 
reduce the jt? cut. To determine to what extent the acceptance can be improved, both the signal 
and background must be studied carefully, because the signal to noise ratio is near unity in this 
mode. Such a study is beyond the scope of this paper. 

We summarize this section as follows. 

• The previously claimed cross section error due to the uncertainty in the acceptance should 
not be taken as a generic statement. The point studied in Ref. [0] has special kinematic 
properties making it a very pessimistic case. 

• Because the error in the acceptance scales inversely with the acceptance, the acceptance un- 
certainty may be minimized by changing the beam energy and the cuts so as to maximize the 
acceptance. One should always try to find the best possible way to increase the acceptance, 
not in order to increase the statistics, but rather to reduce the systematic error. 

• If it is clear that specific distributions and/or cuts are the dominant source of the uncertainty, 
one might benefit from fitting the distribution. Fig. |9] illustrates an example where the shape 
of the acoplanarity angle distribution (not its overall normalization) changes substantially 
with Am. 



the dependence on the VF-boson velocity. The daughter f^-boson in the center of mass frame is substantially 
nonrelativistic at y/s — 400 GeV, so the Vl^-boson velocity is sensitive to the chargino/neutralino mass difference. 
This dependence may be ascertained on an event by event basis. 
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5 Discussion and conclusions 



In this paper we calculated the chargino production cross section including full one-loop quark 
and squark loop corrections. Quark and squark loop corrections are known to induce corrections 
proportional to logm^. This logarithmic correction is seen as a reflection of broken super symmetry 
in the effective theory below the squark mass scale. The correction may be observable in chargino, 
neutralino, and slepton production and decay processes as discussed in Refs. (JTT], [L3|]. In this paper, 
the important corrections in the large mq limit were extracted from the full one-loop calculation, 
and they were compared with the decoupling corrections. We also revisited previous MC studies 
of the measurement of the correction to the fermion-sfermion-chargino coupling Sg = g ei> ^ — 5'f M - 
In Ref. [|TTJ it was stated that a precise determination of the chargino production cross section 
is important, but the uncertainty in the theoretical underlying parameters would be the limiting 
factor in the measurement. In this paper we pointed out the systematic error will not be a problem 
at generic points in MSSM parameter space. This is contrary to their remarks which were based on 
a MC study at a point in parameter space with special kinematic properties. Our study shows that 
experiments at future e + e~ colliders should be sensitive to the squark mass scale if the chargino 
is produced with a large cross section. 

We presented our one-loop calculation in terms of the renormalization scale independent effec- 
tive chargino mixing matrices, U p and V p . They are the matrices which diagonalize the effective 
mass matrix Mc{p 2 ) at momentum p 2 = ?77.|_. When the one-loop amplitude is written in terms 
of U p and V p , a complicated part of the wave function renormalization is absorbed, and the 
remaining part is a simple expression. The sum of the 1PI gauge-x~X + vertex correction and the 
remaining simplified wave function renormalization is scale independent, and decouples in the large 
Tfiq limit. By isolating this scale independent correction we were able to discuss its importance 
separately. 

For sufficiently heavy squarks, it behooves one to introduce the renormalization scale indepen- 
dent effective electron-sneutrino-wino coupling g e -^- All corrections proportional to logm g can 
then be included in the "effective tree level" amplitude, which is obtained by replacing the cou- 
plings and mixing matrices of the tree level amplitude <? e ~^r, U, V, and s-channel gauge coupling 
gi, with the effective ones g ej> ^r, U p , V p , and gf M respectively. The corrections proportional to 
logmq are included in the first three effective parameters, while gf M is rriq independent. 

Since we only include quark and squark loop corrections, only the external chargino lines receive 
wave function renormalization, and we discussed the convenience of the introduction of the on- 
shell effective chargino mixing matrices in that context. We note, however, that our formulation 
of the effective mixing matrices can be easily extended to the gauge-Higgs loops, and to the wave 
function renormalization of other external particles with flavor mixing, such as neutralinos. 

For gaugino-like charginos, the log rriq dependence of g e ^ gives the dominant correction to 
the production amplitude. The amplitude in Higgsino-like chargino production does not receive 
corrections proportional to logm^. Instead it receives finite corrections from the gauge- Higgsino- 
Higgsino vertex correction. The correction is rather small even though Yukawa couplings are 
involved. Numerically we found the correction is of order a few percent. Our numerical calculation 
is in contradiction with previous results given in Ref. They claim large corrections to the 

production cross section of Higgsino-like charginos. We found in some cases order of magnitude 
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differences with their results. Finally, in mixed chargino production, the corrections to the off- 
diagonal elements of the effective chargino mass matrix are important, because U p and V p are 
sensitive to them. The off-diagonal elements receive corrections proportional to logrriq, and they 
also receive decoupling corrections due to squark left-right mixing. Both corrections can be as 
large as 10%. In order to successfully extract the most useful information from the chargino 
measurements, it may be necessary to isolate the light top squark mixing effects, for example 
by measuring the top squark masses and mixing angle through its direct production. Precise 
measurements of the chargino and neutralino spectrum could also give information on top squark 
mixing. 

The validity of various approximations to the one-loop cross section was also studied in this 
paper. A simple approximation which works in a wide region of parameter space makes it easy 
to simulate the effect of the radiative correction in MC studies. We found the 1PI vertex correc- 
tion may be safely neglected, and we further defined approximations to the rest of the one-loop 
amplitude. The one-loop cross section is well described by the effective coupling and mixing matri- 
ces. These parameters encode the leading log(M^) and constant corrections, as well as important 
decoupling corrections. If these decoupling corrections are dropped, we found that the resulting 
approximation can be poor even for relatively heavy squark masses, nig ~ 1.5-^/s. 

Chargino production suffers from IV-boson pair production background at e + e~ colliders. 
Therefore, the detectability of the radiative effect must be studied carefully. Previously stud- 
ies proceeded by choosing a point in MSSM parameter space, and generating the MC signals 
utilizing the cuts that reduce the WW backgrounds while keeping signal events. These cuts were 
determined in a generic situation in Ref. 0. In Sec. 4, we pointed out that the point of param- 
eter space chosen in the MC study of Ref. [|TI]] is not consistent with the assumptions used to 
determine the cuts in Ref. ||. Namely, at the parameter point of Ref. JTTJ there is very little 
phase space in the chargino decay Xi ~^ Wxl- As a result, the $t distribution of the signal 
events is similar to that of the background. However, the cuts to reduce the background were 
chosen under the assumption that the signal would have a higher distribution relative to the 
background. We found this causes the very small acceptance. In this situation small changes in 
the Xi an d Xi masses lead to large variations in the accepted number of events. The expected 
experimental chargino and neutralino mass error therefore leads to a large systematic error in the 
chargino production cross section. We suggest the uncertainty at such a point can be reduced by 
optimizing cuts and beam energy to increase the acceptance. We showed that at generic points 
in parameter space the acceptance is substantially larger, and the systematic errors due to the 
chargino and neutralino mass uncertainties will not pose a serious problem. We stress that efforts 
to optimize cuts to obtain the maximal acceptance greatly reduce the error in the cross section 
both by increasing statistics and reducing systematic errors. This improves the sensitivity of the 
measurement to the loop effects. 

We note the systematic error due to the theoretical underlying parameters may be reduced 
by measuring various kinematical distributions of decay products in ]/> T , 9 acop , etc. Such fitting 
to decay distributions has not been considered in previous studies. Furthermore, the decoupling 
correction is not negligible in the mixed case, and this might introduce an interesting twist in 
future chargino studies. This will be studied elsewhere. We did not present our fits of chargino 
production cross section to MSSM parameters. Notice, however, that fits of MC data to MSSM 
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parameters are sensitive to the specific choice of the theoretical input parameters, beam conditions, 
etc., chosen for the study. The fitted results at one or a few points in parameter space should not 
be interpreted generically. 

The corrections encoded in g ej> ^r, U p and V p are universal. They appear in various production 
and decay processes, and may be important when chargino decay distributions or branching ratios 
are used in a fit. Neutralino pair production receives analogous logm^ corrections. Of course the 
chargino and neutralino corrections are equally important in final states which receive contributions 
from both chargino and neutralino production. 

Previously, information on particles which were not produced directly was ascertained by cal- 
culating the effects of loop corrections in SM processes and comparing the predictions with exper- 
imental data. Unfortunately, superpartners typically give very small corrections in SM processes 
because of their decoupling nature. Once a superpartner is found, the existence of heavier su- 
perpartners with mass M gives rise to interesting non-decoupling effects proportional to log M in 
the production and decay processes of the lighter sparticle. In this paper we studied a chargino 
production process, and compared the log nig correction and the associated decoupling corrections 
in detail. We found the mixing of light third generation squarks also leads important radiative 
corrections. If these two effects can be separated, we could uncover rich information about the 
squark mass spectrum. We stress that a systematic treatment of the loop correction and a detailed 
examination of future experimental prospects are needed to make such a study possible. 
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Appendix A: Tree level interactions 



We list the tree level interactions of charginos, quarks, and squarks. The chargino (x ) mass 
matrix in the gauge eigenbasis, 



is given as 



Vfc = (WE , HH) , V« = (W£, H 2R ) 

-£m = ^mMcijlpjL + 1>iL M Cij*l>7R ' 

M = ( M 2 V2M w cos(3\ 
C \V2M w sm(3 /j ) ' 



(A.l) 



(A.2) 



The mass matrix Mq is diagonalized by two unitary matrices V and U as Md = V*McU\ 
where Md = diag(mj). Note that at the one-loop level My/ in Eq. (A^) is the DR renormalized 
parameter. 
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The chargino-fermion-sfermion couplings are written as follows: 

Ant = -Ki^{a- h ..P L + b- h ..P R ) / 2 + (h.c.) 

-hi W ( a %ij p t + b U p n)h + (h-c-) , (A.3) 

where / = (q,l) and (/i,/2) are SU(2) doublets, and the suffix % of sfermions denote its mass 
eigenstates. Explicit forms of (a, b) in the gauge eigenbasis of sfermions Jl,r are written in terms 
of (U, V), the gauge coupling g 2 and Yukawa couplings yj of fermions / as 

a lu = 92Va , 0,^ = 92^ , 
b] lLl = ~Vh Ui2 , b + hL . = -y fl V i2 , 

b m = > ( A - 4 ) 

where 

= 92m h = g 2 m h 

11 V2M w sm(3 ' /2 V2M w cos/3 ' V ' ' 

The gauge interactions of fermions, charginos, and sfermions are expressed as 

£int = (v%Pl + V&jPr) xJG» 

-h, (vf L P L + v G R P R )fG» 

-ivffif* % (A.6) 

Here G = (7, Z) and 

„£, = (UT G U% , v% = (V*T G V T ) 



Lij — \" j. " Jij > ^ffij — V" - 1 " , 

'V 

1 , _:_2 /1 

VK 



™ f-1 \ ^ z / -l + sin 2 fl 

T =e -1 ' = M -i + sin 2 ^ 



v fL = 9z(T 3fL - Qf sin 9 W ) , = -gzQj sin 6> w , 
/•Z z fz z fz fz n 

V LL = V fL , V RR = V fR , V LR = V RL = , 

v }l = v }R = e Qf , 4 7 = e <?/ % • ( A - 7 ) 

The mixing of left- and right-handed sfermions may not be negligible for third generation 
sfermions. The mass matrices for / = (t, b) are given as follows 

m l ™Ir\ ( h 



Lm -y htR )\ml R m R ){f R J ' 
\ = m~ +m 2 f + m 2 z cos 2(3 (T 3fL - Q f sin 2 W ) , 



m. 



m R = rnj 3R + iTif + ra z Q f cos 2(3 sin 9 W , 

2 _( -m t ( A t + /j, cot (3) fort 
m LR ~ i - mb (A b + //tan/3) for 6 " 1 8j 
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The mass eigenstates /i^ are obtained by diagonalizing the mass matrices. This leads to the field 
rotations 

fx = h cos 9 f - + f R sin 9 f ~ , 

f 2 = -f L dn0 f - + f R coBe f , (A.9) 

where < m? . 6j is the mixing angle. Couplings of /x,2 are easily obtained from those in the 
basis. 

Appendix B: Quark and squark loop functions 

We list the explicit forms of the quark-squark loop functions in the corrected amplitude shown in 
Sec. 2. The results are for quarks and squarks of a given generation. 



The forms of the chargino two-point functions Ey(p ) are [26] 

» dki dkj 



N c 

W > ~ J^^dk^dkj^ 1 ^ ^" b ^" b d h )^^uk^uky 

S S (P 2 ) = ^^(^ki a lj muB ^P 2 '» m «> m 4) + a ^kiKkj m dBo(p 2 , m d , m s J) . (B.l) 

Here _Bo,i are 't Hooft-Veltman functions in the convention of Ref. |T7j], and iV c = 3 is color factor. 

The one-particle-irreducible (1PI) corrections to the xtxjG^ vertices, Fy and F<f, appear in 
Eq. fl6|)ll The corrections have two parts: contributions from (/, /, f' x )-\oops (denoted with /) 
and those from (/, f' x , /y)-loops (denoted with /'), where (/,/') denotes an SU(2) multiplet of 
quarks. Accordingly, the F G 's are decomposed as 

F VL(R) = Y F VL(R) + Y F VL(R) > 
/ / 

F SL(R) = F SL(R) + F SL(R) ■ ( B - 2 ) 
/ / 

The contribution of the (/, /, f' x )-\oops for (/, /') = (d,u) are expressed as 

F vl = TF~2 Y [ b *i V R b *xj FfX + a xi v L a* Xj rrn rrij (Cf 2 x - C(i) 
lt)7r X=l,2 

+m f {-a Xi v R b* Xj rrii C{* - b Xi v L a* Xj rrij C{* 

+a X i v L b* Xj mi (C /x + C/ 2 X ) + b Xi v R a* Xj rrij (C /x + c£*)} 
+m)b Xl v L b X] Cl x ] , (B.3) 



"There we have ignored additional terms proportional to (p3 + Pa)^ since their contribution vanishes in the 
massless electron limit. However, these terms can contribute to other processes such as chargino decays. 
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7°f 
SL 



16tt 2 



X=l,2 



+a Xi v L a* Xj rrij (C{? - C{ 2 + C$? - C(*) 
+m f \a Xi v R b* Xj C{ 2 - a Xi v L b* Xj (C£ x + C(f)j 



Here we abbreviate a Xi — a?, 



(B.4) 

. The Fn^ formula are obtained 



"fxv bxi ~ b f'xv and Vl ( r ) ~ u fMR)- ■ LiiC - 1 R 
from the corresponding F L expressions by replacing a X i <-» b X i and Vl, Vr. Cq and C a g are 
the Passarino and Veltman [27| C functions in the convention of | 28| . The arguments of the C 
function are 



(B.5) 



The function F? x is defined as 

F fX = _1 ^ B fX _ sC fX + ( ^ _ m 2 + m 2 _ s) c fX + (m 2 _ m 2^ _ ^ Q iX _ ^ 



:b.6) 



where B^ x = B (p^,mj, 



in 



fx 



The contributions of (/, f' x , Jy)-loops for (/, /') = (d,u) are expressed as follows 



7 Gf 
VL 



7 Gf 
SL 



16n'- 



2b Xi v X y b* Yj C( 4 



XY 



(B.7) 



XY 



~T^~2 m f a xi vxy b* Yj {C{ 2 XY - C( X XY ) + m, b Xi v XY b* Yj (C£ 3 XY 
ion XY L 

i * (r<PXy rif'XY , r^f'XY ^f'XYs 

+m j a Xi v XY a Yj {Ci 3 - C J 21 + C{ 2 - C[ x ) 



(B.8) 



The are obtained from Eqs. ( |B.7HB.8|) by replacing a Xi <-> b Xi . Here, t>xy = ^xy> anc ^ 



= c o( a /3)(PpPh s > m )> , m p m /' ) • ( B -9) 



The contributions for (/, /') = (u, d) loops can be obtained from the (/, /') = (d, u) expressions 
in Eqs. |||| 



D by replacing a f -, Xi 



b + ~* b~~ 

f'XV f'Xi 



p V f,L -V^fl' alld 



f'G 
} XY- 



Appendix C: Helicity amplitude method 

In the calculation of cross sections it is often useful to directly evaluate the amplitude for helicity 
eigenstates of initial and final state particles, and numerically take the sum of the squared ampli- 
tudes for helicities, instead of taking the trace of the squared amplitude analytically. This method 



is called the helicity amplitude method JT9 



In this appendix we list the spinor bilinears which are relevant in the one-loop amplitude of the 
process e~(pi, hi)e + (p 2 , h 2 ) — > xi (P3-> hsjXj '(P4j ^4) • ^1-4 are the helicities of the corresponding 
particles and take the values ±1/2. We evaluate the spinor bilinears in the center of mass frame. 
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The coordinate space is chosen so that the initial e~ goes along the positive z axis and the final 
Xi goes along the (6 <fi) direction in the polar basis. We give our results in the spherical basis. A 
Lorentz vector in the spherical basis, (A°, A m ) with m — 0, ±, is related to the vector N 1 in the 
Minkowski basis via 



A+ = --j=(A 1 + iA*), 
A° = A 3 , 

A' = 4= (A 1 - %A 2 ) . 

The inner product of two vectors (A, B) is given by 

A ■ B = A B + ^(-l) m+1 A m fi- 



(C.l) 



(C.2) 



The initial massless fermion bilinears H (we ignore the electron mass) and final fermion bilinears 
H are given in the spherical basis as 



Hs 
Hp 



H$ = v h2 (p 2 )Yu hl ( Pl ) = (0,2>/2|A 7 |£<L roiAj ) , 

Uhs(P3)v hi (Pi) = (( E 3 + E 4 ) 2 - (mi + rrij) 2 ) ' 5 Xf o , 
Wht = (-l) h * +1/2 ((E 3 + E 4 ) 2 - K - m 3 ) 2 ) V2 5 Xf0 , 
E 3 + m, yjE 4 - rrij - \J E 3 - m { ^E 4 + 



K = 



E\ EE 



UhsTVh* = 

Uh 3 l^l5Vh 4 = 



rrij 



'x F o, 



E 3 + rrii 
V E 4 + rrij 

_-Qh 4 +l/2 



-E 4 + E 3 + mi + rrij } 5\ f0 , 



' E 4 + rrij 
E 3 + rrii 



(v^)IAfI [ (_i)^+i ^Es-rm ^E 4 + m 3 + yj % + m, yj E 4 - rrij 



dl, F {oy {m - XF)(k 

(C.3) 



We use the abbreviations Xj = h± — h 2 and = h 3 — h 4 . (E, E 3 , E 4 ) are the (e , Xi , Xj) 
energies, respectively. The relevant Wigner rf-functions d^iO) are 



(d\ + d\ d\_\ 

H 1 H 1 H 1 
u o+ u oo u o- 

{d'_ + d\ dl_J 



(0) 



^= sin 6 cos 6 — ^ sin ^ 

1 — COS 1 •_ Q 1+COS0 



(C.4) 
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